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Spontaneous chiral symmetry breaking by hydromagnetic buoyancy 
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Evidence for the parity-breaking nature of the magnetic buoyancy instability in a stably stratified gas is 
reported. In the absence of rotation, no helicity is produced, but the non-helical state is found to be unstable to 
small helical perturbations during the development of the instability. The parity-breaking nature of an instability 
in magnet ohydrodynamics appears to be the first of its kind and has properties similar to those in chiral symmetry 
breaking in biochemistry. Applications to the production of mean fields in galaxy clusters are discussed. 
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The phenomenon of spontaneous breakdown of chiral sym- 
metry, i.e., the bifurcation of an achiral state into two states 
with opposite chirality (handedness), has attracted attention 
in many different fields of physics, even at the macroscopic 
level. A well studied example from hydrodynamics is the 
Taylor-Couette flow with counter-rotating cylinders. In this 
case, linear stability analysis reveals the presence of two op- 
positely helical states (spiral vortices) with identical growth 
rates. Experiments and weakly nonlinear analysis show that, 
depending on the initial conditions, one of these two heli- 
cal states is selected yj]. One might suppose that rotation 
(which is an axial vector) is essential to enable the occurrence 
of such symmetry-degenerate unstable eigenmodes and hence 
the breaking of chiral symmetry. In other words, the conjec- 
ture is that, in the absence of such a vector helical states will 
never show up. Then the question arises naturally whether an- 
other global axial vector can replace rotation in this respect. 
Here the magnetic field is a suggestive candidate. 

However, more puzzlingly, helical flows are possible even 
in non-rotating, non-magnetic setups as demonstrated in 
Ref. 10] for one of the asymmetric square patterns of Benard 
convection. There, the sign of helicity realized in a given cell 
is expected to depend only on the initial condition. Yet, when 
considering that this type of pattern emerges not by a bifurca- 
tion from the trivial (non-convective) solution, but from con- 
vective rolls instead, we can again identify an axial vector in 
the unstable reference state, namely the vorticity of the rolls. 
Note that this vector can, at least with respect to a single roll, 
be considered global. 

As helical flows are known to be crucial in producing 
large-scale magnetic fields by the so-called a effect |3J], the 
strongest interest in such flows is likely to be found in those 
branches of astrophysics, geophysics and planetology which 
deal with the origin of galactic, stellar and planetary magnetic 
fields. However, helicity in astrophysical flows is normally 
not supposed to be caused by a spontaneous parity breaking. 
Indeed, bodies such as the Sun tend to produce helicity with 
opposite signs in their two hemispheres. This is not surpris- 
ing, because the Sun is stratified in the radial direction by its 
own central gravity g and rotates with finite angular velocity 
f2. These two vectors form a pseudoscalar g ■ fl that changes 
sign about the equatorial plane and is, by virtue of the Coriolis 
force, directly related to the kinetic helicity. 



The question now emerges whether finite kinetic and mag- 
netic helicity can also be produced in the absence of ingredi- 
ents that are a priori given as, for example g or fl. In search 
of such a process we have to look for an instability which 
shows a preference for amplifying helical velocity (and pos- 
sibly magnetic) perturbations in comparison with non-helical 
perturbations, but showing of course no preference for one 
specific sign of helicity. One particular example is the mag- 
netic buoyancy instability. It is well known that this instability 
can produce helical magnetic and velocity fields owing to the 
presence of rotation and stratification |4|]. Nevertheless, it has 
been explicitly stated @, Ht] that, in the absence of rotation 
or at the equator, the a effect must vanish because the eigen- 
modes of the instability are then degenerate and have opposite 
helicities. By contrast, recent work [6] delivered indications 
that finite helicity can emerge even at the equator. In this pa- 
per we elaborate on the possibility that helicity may be finite, 
even in the absence of rotation. We present new simulations 
and argue that our findings provide strong evidence for spon- 
taneous chiral symmetry breaking during the nonlinear phase 
of the magnetic buoyancy instability. In addition to our ex- 
ample, there is now also that of the Tayler instability which 
leads to spontaneous chiral symmetry breaking |7|]. In both 
cases there is no rotation and helicity emerges in the absence 
of quantities that are a priori given such as fi and either g or 
at least boundary effects. 

There is quite a number of cases where spontaneous chiral 
symmetry breaking has been discussed. A somewhat hypo- 
thetical example is the production of magnetic helicity during 
electroweak baryogenesis [8] and may be connected with the 
emergence of baryon asymmetry [91]. Chiral symmetry break- 
ing is also known in biochemistry where it refers to the conse- 
quent selection of one of two possible forms of biomolecules 
(mainly sugars and amino acids) that are mirror images of 
each other. This selection may have taken place during the 
emergence of life on Earth lid]. It requires the preferred 
replication of molecules of the same handedness, known as 
autocatalysis ill lL 1 1 211 . Equally important is the effect of mu- 
tual antagonism. This has been identified in the context of 
DNA polymerization, where it characterizes the inability to 
continue polymerization with monomers of opposite handed- 
ness 11311 . Thus, an initial imbalance in chirality is strongly 
amplified and the final selection of one chirality over the other 
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consequently is a result of random fluctuations in the equipar- 
tition of right and left-handed monomers in a prebiotic mix- 
ture 0]. 

The equations governing the agent concentrations in bio- 
chemistry [e.g., Ref. 10] on the one hand and the velocity 
mode amplitudes in hydrodynamics [e.g., [HI] on the other 
show remarkable similarities. Hence, it could be fruitful to 
keep the mentioned biochemical processes in mind when try- 
ing to understand bifurcations into helical states in hydrody- 
namics. In particular, it deserves interest if something like 
'mutual antagonism' can be identified in the nonlinear evolu- 
tion of two competing eigenmodes with opposite helicity. 

We consider a Cartesian slab of ideal gas with a stable strat- 
ification in the z direction which is altered by the presence 
of a narrow horizontal magnetic layer in magnetostatic equi- 
librium. That is, in the layer a fraction of the gas pressure 
is substituted by magnetic pressure. Equilibrium is achieved 
by reducing the temperature correspondingly, but keeping the 
density unchanged, so the magnetic layer is initially not buoy- 
ant. However, the local temperature change in the layer per- 
turbs the thermodynamic equilibrium. We solve the compress- 
ible hydromagnetic equations in a setup described in detail in 
Ref. (a], but focus here on the case without rotation, $7 = 0, 
which shows no a priori preference of positive or negative he- 
licity. 

Our computational domain has an aspect ratio L x : L y : 
L z = 1 : 3 : 1. The ratio of the thickness of the mag- 
netic layer, Hb, to pressure scale height H p and L z is Hb ■ 
Hp : L z = 0.1 : 0.3 : 1. Our fluid and magnetic Prandtl 
numbers are equal to 4. The Lundquist number based on the 
thickness of the magnetic layer is vaqHb /i] = 1000, where 
i>ao = Bo/ y/poHo is the Alfven speed associated with the 
initial magnetic field of strength Bq in the y direction, po is 
the initial density at the bottom of the domain, fiQ the vacuum 
permeability, and r) the magnetic diffusivity. The initial strat- 
ification is a polytrope with index 3, so density scales to tem- 
perature like p ~ T 3 . We use the fully compressible PENCIL 
CODE ll3"(ill for all our calculations and the number of mesh 
points is 64 3 . 

Considering reflectional symmetries of the basic equations 
together with the unperturbed initial state and the boundary 
conditions, we observe that reflections at planes z = const 
clearly change the system already due to its stratification. In 
contrast, reflections at planes y = const leave it unchanged 
while those at planes x — const are equivalent with merely 
changing the sign of the initial field of the magnetic layer. 
Hence there is no essential difference in the behavior of a so- 
lution and a corresponding counterpart obtained by one of the 
two latter reflections. In particular, if there are helical eigen- 
modes of the linearized system they should occur in pairs with 
opposite helicities, but equal growth (or decay) rates. Further, 
if there were stable quasi-stationary helical solutions of the 
nonlinear system they should again occur in such pairs. 

Our results depend decisively on the details of the initial 
velocity perturbations, which possess a slight imbalance of 
random sign in net kinetic helicity. Depending on this sign we 
find exponential amplification of positive or negative kinetic 
and current helicities during the initial stage of the instability. 



An important consequence of this parity breaking is the oc- 
currence of what is known in dynamo theory as the a effect. 
This effect is crucial for the amplification and sustainment of 
a mean magnetic field B. Its evolution is essentially governed 
by the mean electromotive force, £, which is the correlation 
it x 6 of velocity and magnetic field fluctuations, u — U — U 
and b = B — B, respectively, where the overbars indicate hor- 
izontal (xy) averaging. If u lacks reflectional symmetry, and 
in particular, if it exhibits handedness, u x b can have a con- 
stituent proportional to the mean field: 

(u x b)i = aij o Bj - r)ij o Jj, (1) 

where the operator "o" denotes convolution in space, and in 
general also in time, but this will be ignored here. Note that 
and rjij are tensorial kernels, the former a pseudo and the 
latter a true tensor. The symmetric part of ay results in the 
aforementioned a effect, while that of rjij describes what is 
known as turbulent diffusion. Given that we are working with 
horizontal averages, B = B(z,t), we have B z = in the 
absence of an imposed field. Therefore, the evolution of B is 
governed only by 2 x 2 components of both and rjij with 
a xx and a yy being of particular interest. 

In Fourier space (with respect to z) the convolutions in 
Eq. (Q]) turn into multiplications with tensors that depend on z 
wavenumber k. However, due to the intrinsic inhomogeneity 
of our system they depend at the same time also explicitly on 
z. The functions \{z, k) and fjij(z, k) can be directly de- 
termined by the so-called test-field method Ila,ll7ll : for spe- 
cific details see @]. For the purpose of the present paper it 
suffices to consider the Fourier constituents with the smallest 
wavenumber, k = ir/L z and we omit the argument k in what 
follows. 

In this paper we discuss two pairs of runs, A± and B±. For 
each pair, the initial velocity of one run is an exact mirror im- 
age of that of the other one with respect to a vertical plane. 
Consequently, the two runs have opposite initial kinetic helic- 
ity. Depending on its sign, the runs are labeled by + or — . 
In both pairs the initial velocity contains a random pattern of 
horizontal eddies with Mach numbers of the order of 10~ 5 . 
In pair A±, there are additional random vertical motions with 
the same Mach number, whereas in pair B± these are set to 
zero. Thus, the pair A± differs from pair B± in the mag- 
nitude of the initial kinetic helicity, whose normalized value 
(w ■ u) / u rms w Tms is fa 2 x 10~ 5 for the former and, due to 
numerical noise only, ps 4 x 10~ 10 for the latter. In the linear 
phase of the instability the two runs of each pair yield exactly 
the same growth rate. Likewise, the magnetic and kinetic en- 
ergies in the saturated stage are nearly the same. The evolu- 
tion of the normalized kinetic helicity is shown in Fig.Q]for all 
four runs. Here, angular brackets denote the average over the 
full computational domain. Note that parity breaking becomes 
obvious after about one Alfven time, £ao = L v /vaq- In the 
lower panel of Fig. Q]we show the moduli of the kinetic helic- 
ity together with the energy of the components B x z generated 
from the initial B y by the instability. Evidently, the evolution 
of the magnetic energy is nearly the same in all cases. Note 
that, along with the kinetic helicity, the current helicity also 
is generated with the same growth rate, again with positive 
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FIG. 1: (Color online) Upper panel: Evolution of the normalized ki- 
netic helicity for runs A± and B±. Lower panel: Same as above but 
for the absolute value (dotted) overplotted with a solid line to indicate 
the intervals of positive kinetic helicity. Inset: M± — (B x + 13%}, 
normalized by the averaged initial magnetic energy density, Mo = 



B y0 ) for the same set of runs. 



or negative sign, depending on the initial condition. Further- 
more, the final sign of the helicity does not need to agree with 
the initial one. For example, A+ had positive helicity ini- 
tially, but ends up with negative helicity, undergoing several 
sign changes during the kinematic stage. Given that the he- 
licity is a volume integral comprising patches with opposite 
signs of w ■ u, such sign changes can be understood in terms 
of local volume changes of these patches. The sign change 
for the A+ solution at t/t ao = 1.7 could then be explained 
by a gradual shrinking of a patch with positive helicity that 
occupied a much larger fraction at earlier times. We specu- 
late that this scenario might be analogous to the case of chiral 
molecules, where patches of one handedness would shrink in 
the direction of the curvature vector fl8ll . 

In Fig. |2] we plot the dependence of a xx and a yy on z and 
t for the runs A±. Note that their patterns reflect the rise of 
the magnetic layer from its original position at z = 0.3 L z to 
the top of the domain at z = L z . Obviously, within each pair 
of runs the a yy map of the run with negative initial helicity, 
labeled "— ", can be obtained from the one of the run with 
positive initial helicity, labeled "+" to reasonable accuracy by 
a simple sign inversion. 

The component a xx has been found to have the same sign 
throughout the entire domain at any instant of time. Hence it is 
sufficient to consider the evolution of the average (a xx ) which 
is shown in Fig. [3] for all four runs. Its behavior is similar to 
that of (w ■ u) in Fig. Q] which, however, shows a sign change 
at t/tAo ~ 2, whereas the sign of (a xx ) does never change. 
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FIG. 2: (Color online) a xx and a yy as functions of time and height 
for runs A+ (top) and A— (bottom). Values scaled by 10 /vao. 
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FIG. 3: (Color online) (& X x) as a function of time for runs A± and 
B±. 



The phenomenon of spontaneous parity breaking has been 
observed earlier for hydrodynamic instabilities, in particu- 
lar those of the Taylor-Couette flow lfl9ll . For such insta- 
bilities, systems of amplitude equations of weakly nonlin- 
ear, Ginzburg-Landau type have already been described in 
Refs. In all these cases, the helical states, consist- 

ing of traveling waves, emerge due to the presence of rotation 
which is an axial vector. In our case this role is played by the 
magnetic field. Deriving nonlinear amplitude equations for 
the present problem is nontrivial and will be described else- 
where. Here it suffices to say that the amplitude equations in 
Ref. Il5ll show a striking similarity to the equations that de- 
scribe homochirality in the biological context [11]. In particu- 
lar both sets of equations contain terms of mutual antagonism 
which allows one sign of helicity to grow at the expense of the 
other. The analogies between homochirality and the present 
problem extend even further: (a) In both cases one can start 
with an initial state with small excess of one sign of helicity 
or chirality, but the final state can have the opposite sign lfl8ll . 
here, e.g., the set of runs A. (b) In a closed system, which in 
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the biological context implies that the resources available to 
the evolving species are finite, both problems fail to reach a 
statistically stationary state lfl2l. Eoll . 

As eluded to in the beginning, we are not aware of any sim- 
ilar spontaneous parity breaking instability in magnetohydro- 
dynamics. A related example is the selection of one type of 
non-mirror symmetric crystals during their growth in the pres- 
ence of grinding 112 U 12211 . In the case of the hydromagnetic 
buoyancy instability presented here, we can only speculate 
about possible applications. In any case we must be thinking 
of a stably stratified layer with negligible rotation. A prime 
example fitting this description is clusters of galaxies. Indeed, 
magnetic and thermal buoyancy effects are invoked to explain 
what looks like cold bubbles in such clusters [23]. 

However, the magnetic buoyancy instability might be just 
one example showing the parity-breaking property described 
here. The magneto-rotational instability (MRI) could be an- 
other one. In that case there is rotation, but in the absence 
of stratification and current density there would still be no 
pseudo-scalar constructible and hence no net helicity, unless 
the MRI were capable of spontaneous parity breaking. 

Large-scale magnetic fields have been found in simulations 
without net helicity including cases with shear, in particular 
the MRI 12411 and the magnetic buoyancy instability with a 
sinusoidal shear profile 02511 . Magnetic field generation oc- 
curs also in the case of a linear shear flow with isotropic non- 
helically forced turbulence Although it is not clear whether in 
this case an instability (like the magnetic buoyancy instabil- 
ity) is really required for parity breaking, it should be noted 
that there is generally a production of mean vorticity rf26il27tl . 



whose possible importance has been discussed in connection 
with dynamos from linear shear flows B28I1 . Such a vorticity 
production would indeed be the result of a (hydrodynamic) 
instability, supporting the idea that an instability might be 
needed for parity breaking. An alternative explanation for the 
origin of large-scale fields in this setup is the incoherent a— 
shear dynamo ll29ll . where kinetic helicity is expected to fluc- 
tuate about zero. Revisiting this assertion in the light of the 
present results could be a worthwhile effort. 

In conclusion, we have found a mechanism that is able to 
amplify magnetic and velocity fields with either sign of he- 
licity, but the selection of one of them occurs during the non- 
linear stage where it gives rise to an a effect under otherwise 
parity-invariant conditions. This mechanism might be appli- 
cable to clusters of galaxies. Moreover, it is very likely that 
the magnetic buoyancy instability is just a first example in a 
class of spontaneous parity-breaking instabilities in magneto- 
hydrodynamics. 
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